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Abstract 

In gravitational Higgs mechanism graviton components acquire mass via 
spontaneous diffeomorphism breaking by scalar vacuum expectation values. 
We point out that in the massless limit the resulting theory is not Einstein- 
Hilbert gravity (EHG) but constrained gravity (CG). Consequently, massive 
solutions in the massless limit must be compared to those in CG (as opposed 
to EHG). We discuss spherically symmetric solutions in this context. The 
Schwarzschild solution in EHG can be coordinate-transformed such that it 
is also a solution in CG. The resulting solutions are non-perturbative in the 
asymptotic regime, and are reproduced in the massless limit of asymptotic 
massive solutions, hence no van Dam-Veltman-Zakharov discontinuity. We 
point out that higher curvature terms must be included to obtain non-singular 
spherically symmetric massive solutions and discuss a suitable framework. 



1 Introduction and Summary 



In gravitational Higgs mechanism pQ graviton components acquire mass via sponta- 
neous breaking of diffeomorphisms by scalar vacuum expectation values^. Recently, 
gravitational Higgs mechanism was revisited in the context of obtaining massive 
gravity directly in four dimensions [TBI [T71 \TE\ [TP"] . A general Lorentz invariant mass 
term for the graviton Hmn is of the form 



where (3 is a dimensionless parameter. For (3^1 the trace component h^, which is 
a ghost, is a propagating degree of freedom, while it decouples for the Pauli-Fierz 
mass term with (3 = 1. Massive gravity in Minkowski space with (3 = 1 can be 
obtained via gravitational Higgs mechanism by including higher derivative terms in 
the scalar sector and appropriately tuning the cosmological constant [T9j . 

The framework of [19] provides a ghost-free, non-linear and fully co variant def- 
inition of massive gravity via gravitational Higgs mechanism with spontaneously 
(as opposed to explicitly) broken diffeomorphisms. In this paper we discuss the 
massless limit of gravitational Higgs mechanism within this framework. Our key 
observation is that in the massless limit the resulting theory is not Einstein-Hilbert 
gravity (EHG) but constrained gravity (CG). Consequently, massive solutions in the 
massless limit must be compared to those in CG (as opposed to EHG). For exam- 
ple, the usual Schwarzschild solution in EHG does not satisfy the constraint in CG, 
so we cannot compare massive solutions in the massless limit to the Schwarzschild 
solution in EHG. On the other hand, the Schwarzschild solution in EHG can be 
coordinate-transformed such that it is also a solution in CG. The resulting solutions 
are non-perturbative in the asymptotic regime, and are reproduced in the massless 
limit of asymptotic massive solutions. Therefore, there is no van Dam-Veltman- 
Zakharov (vDVZ) discontinuity [2UJ [2T] in the massless limit of gravitational Higgs 
mechanism in the spirit of [22p . 




We also study massive solutions in the limit corresponding to the near-horizon 
regime in massless solutions. We point out that in the lowest-curvature approxima- 
tion (Einstein-Hilbert term) massive solutions have a naked singularity. This implies 
that the lowest-curvature approximation breaks down, and higher curvature terms 
must be included to obtain non-singular solutions. Adding higher curvature terms 
does not modify the no vDVZ conclusion in the asymptotic regime. We discuss a 
framework for constructing non-singular solutions for all values of the graviton mass 
in the context of gravitational Higgs mechanism. 

iFor earlier and subsequent related works, see, e.g., [HE1S1IS1IS1IZIISIISI HDl EH H21 H3I - For 
a recent review of massive gravity in the context of infinite volume extra dimensions, see |14j and 
references therein. For a recent review of spontaneous breaking of diffeomorphism symmetry in the 
context of Lorentz violating Chern-Simons modification of gravity, see [15] and references therein. 

2 In the context of brane induced gravity no vDVZ discontinuity was argued in [531 HH 
by studying non-perturbative solutions. 
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2 Gravitational Higgs Mechanism 



In this section we review the construction of [19J, and discuss the massless limit of 
gravitational Higgs mechanism. We consider gravity in D dimensions coupled to 
scalar fields <fi A , A = 0, . . . , D — 1. We will take the metric Zab for the global scalar 
indices A,B,... to be the Minkowski metric: Zab — Vab- 

Coordinate-dependent scalar VEVs break diffeomorphisms spontaneously. This 
results in massive gravity. To obtain a unitary theory (i.e., the Pauli-Fierz mass 
term for the graviton in the linearized theory), we must include higher derivative 
terms in the scalar sector. In fact, adding four-derivative terms suffices, but for now 
we will keep our discussion here more general. 

Thus, consider the following action: 

S Y = Mp~ 2 J d D xV^G [R - V{Y)} , (2) 

where a priori V(Y) is a generic function of Y, and 

Ymn = Z AB V M 4> A V N 4> B , (3) 
Y = Y MN G MN . (4) 

The equations of motion read: 

V A/ (V'(F)V M A ) = , (5) 
Run - \GmnR = V\Y)Y MN - ^G MN V{Y) , (6) 

where prime denotes derivative w.r.t. Y. Multiplying (jHJ) by Zab^s<P B and con- 
tracting indices, we can rewrite the scalar equations of motion as follows: 



d 



M 



V^GV'(Y)G MN Y NS ] - l^GV'(Y)G MN d s Y MN = . (7) 



Since the theory possesses full diffeomorphism symmetry, (J7J) and ([6]) are not all 
independent but linearly related due to Bianchi identities. Thus, multiplying dSJ) by 
y—G, differentiating w.r.t. V N and contracting indices we arrive at (171). 
We are interested in finding solutions of the form: 

^ = m 5 A M x M , (8) 
Gmn = Vmn ■ (9) 

These exist for a class of "potentials" V(Y) such that the following equation 

Y V\Y) = | V(X) (10) 
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has non-trivial solutions^. Let us denote such a solution by Y*. Then we have 



m 2 = Y*/D . (11) 

Note that (jTOl) is not a tuning condition on couplings in V. Rather, it fixes the 
scalar VEVs (whose slopes are given by m) in terms of various couplings in V . 

The equations of motion §5$) and (JH]) are invariant under the full diffeomorphisms 
of the theory (G MN = i] MN + h MN ): 

5 V A = Wu^e 1 = m 5 A M e 1 , (12) 
Sh MN = V m Cjv + Vat^m • (13) 

The scalar fluctuations ip A around the background ([8]) can be gauged away (set to 
zero) using the diffeomorphisms ffl2|) . Once we gauge away the scalars, diffeomor- 
phisms can no longer be used to gauge away any of the graviton components h MN . 
After setting tp A = 0, in the linearized approximation we have {h = hj^): 

d N h MN -(3d M h = , (14) 
1 Ml 

Rmn — -^GmnR = —y- (PvMNh — Iimn) , (15) 

where 

q = - -im^"^ (16) 

Ml = V(Y) . (17) 

In particular, for a special class of potentials with 

4Y # V'(Y*) = -D 2 V(Y*) , (18) 

we have (3 = 1, and the correct tensor structure for massive gravity without non- 
unitary states. In this case the explicit form of the equations of motion is given 
by: 

h = 0, (19) 
d N h MN = , (20) 
d s d s h MN = M 2 h h MN , (21) 

In particular, the ghost state h decouples, and we have massive gravity with (D + 
1)(D — 2)/2 propagating degrees of freedom. 

Thus, we can get the Pauli-Fierz combination of the mass terms for the graviton 
if we tune one combination of couplings. In fact, this tuning is nothing but tuning of 



3 This condition is not particularly constraining. Such solutions exist even for linear potentials 
V = A + Y. In this case the restriction is that the cosmological constant be negative. 
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the cosmological constant - indeed, ([TBI) relates the cosmological constant to higher 
derivative couplings. 

Thus, consider a simple example: 

V = A + Y + \Y 2 . (22) 

The first term is the cosmological constant, the second term is the kinetic term for 
the scalars (which can always be normalized such that the corresponding coefficient 
is 1), and the third term is a four- derivative term. We then have: 

y* = — r A- 1 , (23) 

2(D + 2) V ; 

which relates the mass parameter m to the higher derivative coupling A: 

m2 = y */ D =-2pT2) A "' (24) 
The graviton mass in this case is given by: 

M * = -JdTW A " ■ (25) 

Note that we must have A < 0. Moreover, the cosmological constant A must be 
tuned as follows: 

A = D2 ; 4D ~ S A- , (26) 

4(£> + 2) 2 ' v ; 

which implies that the cosmological constant must be negative. 



2.1 Constrained Gravity as the Massless Limit 

Here we would like to begin addressing the question of what happens in the limit 
where the graviton mass goes to zero. A naive expectation that one recovers 
Einstein-Hilbert gravity in this limit does not hold. Instead, we get constrained 
gravity. 

To see what happens in the massless limit, for illustrative purposes, let us con- 
sider the example with the quartic potential given by (l22l . (The analysis in the 
general case completely parallels our discussion here, so we will not repeat it.) To 
obtain the massless graviton limit, we must take A — > — oo. In this limit both Y* 
and A go to zero, and so does m 2 . Let us gauge fix the scalars to zero via (|SJ). We 
then have: 

Y MN = m 2 E MN . (27) 

Here Emn is the Minkowski metric in the coordinate frame given by x M . It coincides 
with the flat Minkowski metric tjmn if x are Minkowski coordinates. However, in 
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the general case the metric Emn is n °t flat. For instance, in spherical coordinates 
we have 

E MN dx M dx N = -dt 2 + dr 2 + r 2 lab dx a dx h , (28) 

where ^ a b is a metric on the unit sphere S d_1 , d = D — 1. 

In the massless limit we then have the following equations of motion: 







M 



V^GQG MN E NS \ - -^GQG MN d s E MN = , (29) 
Rmn — -GmnR — , (30) 



where 



Q = 1 - ^ G MN E MN . (31) 

Note that fl29|) and fl30|) are not linearized but exact. Eqn. (130)) is just the Einstein 
equation of motion without the cosmological constant. However, gravity in the 
massless limit is not Einstein-Hilbert gravity as we have a constraint (|2"9~|) . 

The fact that we obtain constrained gravity in the massless limit is important. 
In particular, if we take, say, a spherically symmetric solution in massive grav- 
ity and consider the massless limit, we should not expect it to coincide with the 
Schwarzschild solution of Einstein-Hilbert gravity. That would be like comparing 
apples to oranges. Instead, it should coincide with the corresponding spherically 
symmetric solution in constrained gravity. The latter will look very different from 
the Schwarzschild solution of Einstein-Hilbert gravity. One way to construct so- 
lutions in constrained gravity is to start with known solutions in Einstein-Hilbert 
gravity and coordinate-transform them to satisfy the constraint (this is similar to 
[2"T]). In other words, to compare apples to apples, we can, for example, compare the 
massless limit of a spherically symmetric solution in massive gravity to an appro- 
priately coordinate-transformed Schwarzschild solution of Einstein-Hilbert gravity 
such that it satisfies the constraint. 



3 Spherically Symmetric Solutions 

Before we discuss explicit solutions, let us cover some generalities. We will be 
interested in spherically symmetric solutions to massive and massless equations of 
motion (A, B, C are functions of r only): 

ds 2 = -A 2 dt 2 + B 2 dr 2 + C 2 -f ab dx a dx b . (32) 

In the massive case, the Einstein equations read (for definiteness, we focus on the 
quartic potential (|22l) ): 

Rmn = m 2 QE MN + -^—^ [A - XY 2 ] G MN , (33) 
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where Q was defined in (13TI) . Y is given by 

Y = m 2 [A- 2 + B- 2 + (D — 2)r 2 C- 2 ] 



(34) 



and we gave A, A and m 2 in the previous section. 

The non- vanishing components of Rmn are given by (prime denotes derivative 
w.r.t. r): 



Roo 



A B~ 



A AB 
A" A'B 1 
~A ~ ~AB 

Rab = —labR* , 



A " A ' B ' {D -2) A ' C ' 



(D-2) 



AC 
C" B'C 
~~C ~ ~BC 



R„ 



C" 



c 



a 



A^_B^_ 
~A~~B 



The scalar equations of motion (£7§ reduce to a single equation: 
d r [AB- l C D - 2 Q] — (D — 2)rABC D - 4 Q = 



(35) 

(36) 
(37) 

(D-3). (38) 
(39) 



We will now discuss solutions to these equations. In D = 4 we have some simplifi- 
cations, so we will focus on D = 4 for the remainder of this paper. 

3.1 Four-dimensional Massless Solutions 

Let us start by analyzing the above equations in the massless case (m 2 = 0) in 
D = 4. We have the following equations 



(40) 
(41) 

(42) 



A" 
~A 


A'B' 
AB 


A'C 
+ 2 AC 




o, 


A" 


A'B' 




"C" 




B'C 


~A 


AB 


+ 2 


~c 




BC 


C" 


<?) 


2 


c 


'A' 


B' 


~C 


+ 


~c 


~A 


~ ~B 



B 2 C- 



plus the constraint 
where 



d r [AB _1 C 2 Q] - 2rABQ = 



Q = 1 - - \A~ 2 + B- 2 + 2r 2 C~ 2 \ 
6 



(43) 
(44) 



If it were not for the constraint, we could simply take the Schwarzschild solution: 

A = B- 1 = ^l- r f, (45) 
C = r . (46) 
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However, this solution does not satisfy the constraint. 

Actually, there is a systematic way of fining solutions that satisfy the constraint 
by transforming known solutions that satisfy unconstrained Einstein's equations. 
Thus, we start from a known unconstrained solution given by A, B, C, and transform 
the radial coordinate r — > f(r). The resulting warp factors are given by 



A{r) = A(f(r)) , 
B(r) = B(f(r))f(r) 
C(r) = C(/(r)) , 



(47) 
(48) 
(49) 



and they still satisfy the equations of motion. This is because the massless equations 
of motion possess full reparametrization invariance. The constraint then produces 
a second order differential equation for the function f(r). Thus, starting with the 
Schwarzschild solution, we can obtain solutions satisfying the constraint by setting 
f(r) = C(r) in the above expressions, which gives a differential equation for C. We 
have: 



y/l - n/C ' 
and the differential equation for C reads: 

d r [A 2 C 2 Q/C] - 2rQC = . 



(50) 
(51) 

(52) 



While (1521) is highly non-linear, we can solve it in two regimes: near the horizon 
(C — > r*), and asymptotically (r ^> r*). 



3.1.1 Asymptotic Behavior 

In the asymptotic regime we have r ^> r*. We can find a solution via the following 
Ansatz: 



C = r 



r J r V r 



(53) 



where a and f3 are numerical coefficients to be determined. This indeed solves fl52l) 

r* 



A = l- — + 

2r \ r 



v 3 
T*\ 2 



B 



C 



39 V r J 2r \r J ' 
39 V r / r \ r 




(54) 
(55) 
(56) 
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and (3 is an integration constant. This is because we started with the Schwarzschild 
solution and transformed it via r — > C(r). The constraint (I52p is a second order 
differential equation for C, whose solution contains two integration constants. How- 
ever, because we drop subleading terms, the resulting equation actually is effectively 
only a first order equation for C, so we have one integration constant (and the sec- 
ond integration constant controls the subleading terms). It simply parameterizes 
the Schwarzschild solution in the transformed coordinate frame. 

Note that the above asymptotic solution is non-perturbative in the following 
sense. If we write 

C = r(l + c), (57) 

in the asymptotic regime c is small, so we can try to solve (1521) in the linearized 
approximation. However, one finds no solution when keeping only linear terms, so 
one must also keep quadratic terms in c. 

3.1.2 Near-horizon Behavior 

In the near- horizon regime we have C — > r*. In fact, in this case we have r — > 0, 
and: 

C = r* + 7 r + Oir/r,) 2 , (58) 
A = vV/r* [1 + 0(r/n)] , (59) 
B = VV.A* t 1 + 0(r/n)] . (60) 

These warp factors satisfy (1521) for any positive integration constant 7. The near- 
horizon metric reads: 

ds 2 = -^dt 2 + ^—dr 2 + (r* + ~ir) 2 lab dx a dx h . (61) 

This is the metric we ought to compare the massless limit of massive solutions to, 
and not to the Schwarzschild metric, in the near-horizon regime. 



3.2 Four-dimensional Massive Solutions 

Let us now go back to the massive equations of motion. In four dimensions we have: 



A 2 B~ 2 



A" 
~A 



'A" 
_~A 
A'B' 



A'B' A'C 
+ 2- 



AB 



+ 2 



AB 
C" 
~C 



C" (C 



C I B- 2 {^—+\^- \ 



C 



c 



AC 
B'C 
BC 

C 



m 



m 



A 
B 2 



24 



6 



" X 2 " 




X 


~ 24. 




1 

6 



c 



A!_ 
T 



B' 
~B 



m 















V-- 


-r 2 


fl-- 


} 




24 _ 




6 _ 





(62) 
(63) 



(64) 
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plus the scalar equations of motion: 

d r [AB^C^Q - X)] - 2rAB(Q - X) = , (65) 

where 

X = A~ 2 + B- 2 + 2r 2 C~ 2 . (66) 

Just as in the previous subsection, here we are interested in studying solutions 
to these equations in the asymptotic and near-horizon regimes. We will discuss 
perturbative asymptotic solutions first. 

3.2.1 Perturbative Asymptotic Solutions 

Here our goal is to discuss the vDVZ discontinuity in perturbative asymptotic solu- 
tions. Let 

A = 1 + a , (67) 
B = 1 + b , (68) 
C = r(l + c). (69) 

Here we assume that a, b, c go to zero asymptotically, and in the equations of motion 
we keep only linear terms in a, b, c. As we will see in a moment, this approximation 
breaks down for small graviton mass, hence the vDVZ discontinuity. 

The above four equations of motion in the linearized approximation read: 

9 m 2 

a" + -a' = — [3a + b + 2c} , (70) 

T O 

a" + 2c" + -d --b' = ^-[a + 3b + 2c] , (71) 
r r 3 

4 1 2 m 2 

c" + -c' + -[a'-b']--[b-c] = — [a + b + 4c] , (72) 

a' + 2c' - - [b - c] = . (73) 
r 

The solution is given by: 

a" + - a ' = fx 2 a , (74) 
r 

b = -4~> ( 75 ) 
fi z r 

c=-i(a + 6), (76) 



where 



2 

» =-z- ( 77 ) 
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is the mass squared of the graviton. Note that we have 



a + b + 2c = , (78) 

which is simply a manifestation of the fact that the trace of the graviton is not a 
physical degree of freedom. 
Thus, we have: 

a = -e"^ , (79) 

r 

6=-^[l+/,r]e--, (80) 

c=-^[l + /ir + / ,V]e-^, (81) 

where ( is an integration constant. The only way to have a smooth massless limit 
would be to take fi — » and ( — > with ICI/yU 2 = r \ fixed. In this case we would 
have vanishing a in the massless limit, i.e., only b and c would be non- vanishing 

a = , (82) 
b = -2c = -| . (83) 

However, the corresponding metric is equivalent to a coordinate-transformed flat 
metric (in spherical coordinates). 

Above we assumed that A, B and C/r all asymptotically go to 1. In fact, there 
are no asymptotic solutions other than those where A, B and C/r all go to 1. To 
see this, let us assume that A, B and C/r asymptotically go to some constant values 
A, B and C. Eqn. (|65|) then implies that B = C, and the equations ([621 and (|63|) 
imply that A = B. Noting that asymptotically we then have X —>■ 4/ A 2 , Eqn. (l62j) 
implies that A = 1. 

Thus, as we see, perturbative asymptotic massive solutions do not have a smooth 
massless limit (except for trivial solutions equivalent to coordinate-transformed flat 
solutions). This is the well-known vDVZ discontinuity [20],[2T]. However, this dis- 
continuity is an artifact of the perturbative approximation, which breaks down at 
r ~ 7*1. Note that |£| is expected to be of order of the Schwarzschild radius r*, so we 
have 7*i ~ (r*//i 2 ) 1//3 . This scale goes to infinity when /i goes to zero, so one must 
consider non-perturbative solutions [22]. 

3.2.2 Non-perturbative Asymptotic Solutions 

The correct way to think about the massless limit is therefore to consider non- 
perturbative massive solutions. In the massless limit they smoothly go to the asymp- 
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totic massless solutions we discussed in Subsection 3.1.2. Thus, we have: 




C = r 



1 + vl(7) l+/3 7 + °(7) f + ^ 2 ^) 



(84) 
(85) 
(86) 



where (3 is an integration constant. Note that the expansion in /i 2 is valid at distance 
scales r <C l//i- As /i — > 0, we have a smooth massless limit for all r. 

3.2.3 Singular Solutions 

Let us finally comment on what happens in massive solutions in the analog of the 
near-horizon regime in the massless case. Recall that in the massless case at the 
horizon, which in the coordinate-transformed Schwarzschild solution is located at 
r — > 0, we had A — > 0, B — > oo, C — ► r*, with vanishing scalar curvature. In the 
massive case the situation is different. The scalar curvature is given by: 

R = m 2 (X-A). (87) 

If A goes to zero, then X diverges, and so does R, which implies that if A goes to 
zero, we have a naked singularity instead of a horizon. 

Here we can try to look for solutions where A goes to a constant instead. In this 
case X goes to a constant as well, so the scalar curvature is finite. However, here 
we argue that such solutions cannot reproduce the massless near-horizon behavior 
in the massless limit. 

Thus, let us assume that in the massive case A goes to a constant A at the 
horizon. To have a smooth massless limit we must assume that Aq goes to zero 
in the massless limit, and Xq, the value of X at the horizon, then goes to infinity 
(which is the behavior we have in the massless case). For sufficiently small graviton 
mass we can then assume that X ^ 6. 

Next, ([55]) can be written as: 

A' B' 2C X' B 2 

i-B + --(rx = 2r c^ (88) 

In the near-horizon regime we are assuming that C also goes to a constant Co. We 
then have the following leading behavior for B: 

B -^, C ( l V ( 89 ) 
2a/ W - r o) 

where r > is the location of the horizon. Recall that in the massless case the 
horizon is located at r = 0. If we assume that ro = 0, then we have B divergent as 1/r 
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instead of l/\/r (which is the near-horizon behavior in the massless case). Moreover, 
according to (1891) we cannot assume that r goes to zero in the massless limit. 
Therefore, solutions where A goes to a constant at the horizon cannot reproduce 
the massless solutions in the massless limit. 

4 Non-singular Massive Solutions 

As we saw in the previous section, there is a smooth massless limit in the asymp- 
totic regime in the massive spherically symmetric solutions. On the other hand, 
massive spherically symmetric solutions are sing ulaJJ This simply means that the 
approximation we are using breaks down. In particular, appearance of a naked 
singularity implies that higher curvature terms become important and should be 
included. Conceptually, this is similar to the appearance of naked singularities in 
higher codimension brane solutions in Einstein-Hilbert gravity [29] . Once we include 
higher curvature terms, we expect to find smooth solutions, in particular, by using 
the smoothing out procedure of [30] . In the context of higher codimension branes 
explicit non-singular solutions were indeed found in [31j. 

Similarly, we expect that non-singular spherically symmetric solutions also exist 
in the context of massive gravity once we add higher curvature terms. The latter 
generally make the problem technically challenging. However, there are special 
higher curvature combinations, in particular, (lower-dimensional) Euler invariants, 
one can add that are sufficient to get rid of singularities, yet are simple enough to 
keep the problem computationally tractable. 

Thus, we can add a Gauss-Bonnet combination of quadratic in curvature terms, 
whose simplifying property is that all equations of motion still remain second or- 
der differential equations and no third or forth derivatives appear. However, the 
Gauss-Bonnet combination is topological in D = 4, so for it to contribute we must 
consider D > 4. In fact, in D = 5 completely smooth spherically symmetric ana- 
logues of the Schwarzschild solution were explicitly constructed in [32] for massless 
Einstein-Hilbert-Gauss-Bonnet gravity (EHGBG) in D — 5 in a different context. 
We therefore expect non-singular spherically symmetric solutions in massive (via 
gravitational Higgs mechanism) EHGBG in D = 5 to exist, and perhaps even be 
computationally tractable. Note that adding higher curvature terms only affects 
high curvature regions, so it does not modify the no vDVZ conclusion in the asymp- 
totic regime, while it should allow for a smooth massless limit at all distance scales. 
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